Antarctic Radio Frequency Albedo and Implications for Cosmic Ray Reconstruction 
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Abstract 



00 22 From an elevation of ~38 km, the balloon-borne ANtarctic Impulsive Transient Antenna (ANITA) is designed to detect the up- 

23 coming radio frequency (RF) signal resulting from a sub-surface neutrino-nucleon collision. Although no neutrinos have been 

I p4 discovered thus far, ANITA is nevertheless the only experiment to self-trigger on radio frequency emissions from cosmic-ray 

25 induced atmospheric air showers. In the majority of those cases, down-coming RF signals are observed via their reflection from 

I"! 26 the Antarctic ice sheet and back up to the ANITA interferometer. Estimating the energy scale of the incident cosmic rays therefore 
requires an estimate of the fractional power reflected at the air-ice interface. Similarly, inferring the energy of neutrinos interacting 

Q_28 in-ice from observations of the upwards-directed signal refracting out to ANITA also requires consideration of signal coherence 

I 29 across the interface. By comparing the direct Solar RF signal intensity measured with ANITA to the surface-reflected Solar signal 
intensity, as a function of incident elevation angle relative to the surface 6*,-, we estimate the power reflection coefficients 'RiOi). 
We find general consistency between our average measurements and the values of HiOi) expected from the Fresnel equationsHI, 

32 separately for horizontal- vs. vertical-polarizations ("HPol" and "VPol", respectively). 

33 Keywords: radioglaciology, interferometry, cosmic rays 

> 

p^34 1. Introduction 

^^35 1.1. The ANITA Experiment 

* 36 Initiated in 2003, the Antarctic Impulsive Transient Antenna (ANITA) is a balloon-borne antenna array primarily designed to 

(■-^ 37 detect radio wave pulses caused by neutrino collisions with ice[2J. The basic instrument consists of a suite of 40 quad-ridged 

38 horn antennas, optimized over the frequency range 200-1200 MHz, with separate outputs for vertically vs. horizontally-polarized 

^-H 39 incident radio frequency signals, mounted to a high-altitude balloon. From an elevation of ~38 km, the payload observes the 

^ 40 Antarctic continent in a circumpolar trajectory. Two one-month long missions (ANITA-I; Dec. 2006-Jan. 2007 and ANITA-II; 

• '~j 41 Dec. 2008-Jan. 2009) have yielded world's-best limits to the flux of Ultra-High Energy ("UHE", corresponding to energies in 

rN 42 excess of 10'^ eV) neutrinos in the energy range to which ANITA is sensitive p],l5l|. Nevertheless, model-dependent calculations 

^ 43 of surface effects show considerable variation in the fraction of signal power emerging from within the dense ice target medium 

44 and across the interface to the receiver; this is also true for experiments seeking to measure neutrino interactions within the lunar 

45 regolith|l6l |7] |8] 121. An interferometric analysis of the ANITA-I data sample provided a statistically large (16 events) sample 

46 of radio frequency signals attributed to the geomagnetic + Askaryan radiation associated with cosmic-ray induced extensive air 

47 showers lfT2ll (EAS); those events marked the first self-triggered detection of EAS via their radio emissions. 

48 The analysis described below is based on ANITA-II data. After launching from McMurdo Station in December, 2008, ANITA- 

49 II was aloft for a period of 31 days with a typical instantaneous duty cycle exceeding 95%. Although the external experimental 

50 conditions (Uvetime and flight path) were more favorable for ANITA-II than for ANITA-I, the sample of detected extensive air- 

51 shower (EAS) events was somewhat smaller for ANITA-II, primarily owing to a modification to the trigger that favored neutrino 

52 detection, but, as realized ex post facto, dis-favored EAS detectionL10|. The ANITA-III launch, scheduled for December, 2013, 

53 is expected to result in an improvement by a factor of ~100 in the number of detected air shower events relative to ANITA-I. 

54 Correspondingly, we have sought to determine the surface reflectivity in advance of that upcoming ANITA-III mission. 
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1.2. ANITA Data Acquisition and Signal Reception Chain 



The ANITA detector is designed to efficiently collect radio signal over a 1 GHz bandwidth, with electronics and trigger config- 
ured to produce the minimum possible energy threshold. The front-end Seavey quad-ridge antennas provide the first element in 
the ANITA radio wave signal processing chain. These horn antennas have separate VPol vs. HPol signal polarization readouts, 
and each has a field-of-view of approximately 25° half-width-half-maximum (FWHM) in both azimuth angle (0) and zenith angle 
(9). Following the Seavey antennas, high- and low-pass filtering restricts the system bandpass to 200-1200 MHz, after which 
signals are split into a 'trigger' path consisting of a heirarchy of trigger conditions ("levels") and a 'digitization' path. If all levels 
of the three-tiered, heirarchical trigger are satisfied, the digitized signals, consisting of 256 voltages sampled for each antenna at 
2.6 GSa/s, are then stored to disk. During typical ANITA-II data-taking conditions, data were written at a rate of approximately 
10 Hz. The overwhelming majority of these triggers are due to upward excursions of the thermal radio photosphere. 

1.3. ANITA interferometry 

Beginning with the eighty 100-ns duration voltage waveforms which are taken during an ANITA event trigger (40 antennas X 2 
polarizations per antenna), familiar radio interferometry techniques are used to identify localized radio sources. To determine the 
direction of an external impulsive source, the full solid angle is first divided into an Ng x A^,^ grid (typically, 180 x 360) in elevation 
and azimuth, respectively. Next, an inter-channel cross-correlation sum is calculated for each point in the grid, by shifting cross- 
correlated waveforms by the relative propagation time delays that would be incurred, channel-by-channel, for a putative source at 
that (6,<p) grid point. Formally, the cross-correlation C, for two channels ; and j can be written as the dot product of the voltage time 
series for channel i V,(f) and the voltage time series for channel j Vj{t), shifted by the calculated time delay 6t^j(6, (fi) corresponding 
to the putative source S at {6,(p), for the N samples in the dot product (equal to the number of samples in the waveform minus the 
number of time samples corresponding to the time shift dtij). To minimize sensitivity to channel-dependent gain variations, this 
cross-correlation coefficient is normalized relative to the rms voltages in the two channels cr(y,) and o-{Vj): 



In the limit of a strong source which dominates the 100-ns waveform, the cross-correlation rises linearly with signal voltage , 
since cr{V) ~ . In the limit of a very weak source, cr(y) is largely independent of and the cross-correlation increases 
quadratically with signal voltage. The high signal-to-noise (SNR) events that comprise the ANITA EAS sample correspond to the 
former case. The Solar signal described below corresponds to the latter case, and therefore, in the case of signals from the Sun, 
the cross-correlation coefficient scales with the radio power received. The total cross-correlation C is calculated by summing over 
all independent (i,j) pairs. The polar and azimuthal angles of a putative source are then assigned to that grid pixel giving the 
highest total cross-correlation C. More extensive details on the construction of the ANITA interferometric map can be found in 
the UteraturelfTTII. 

An obvious design goal of ANITA is to localize sources with the best-possible pointing resolution. Unlike source vertexing 
algorithms using hit-time assignments from threshold-crossings, interferometry is particularly suited for continuous-wave (CW) 
sources. As illustrated by Monte Carlo simulations (Fig. [TJ, the interferometric reconstruction resolution improves, roughly 
linearly, with the frequency of the source. Alternately, the angular size of the interferometric image provides information on the 
source frequency content, with wider images implying power spectra weighted preferentially to lower frequencies. 

A discone transmitter, buried approximately 100 m into the ice sheet at a remote site near Taylor Dome, has been used to assess 
the inteferometric source reconstruction resolution for signals with broad frequency content. For that transmitter, we obtain values 
of 0.25 degrees in elevation and 0.56 degrees in azimuth for the ANITA-II flight. For comparison, the intrinsic Solar radio source 
size is of order 1 degree in the sky. 

1.4. ANITA Radio Air Shower Observations 

A maximum likelihood analysis of the 16 ANITA-I EAS detections implied a mean energy of 1.5±0.4(stat);^Q3(sys)xlO'^ 
eV for the cosmic ray primary! 12J. with the highest primary estimated at an energy approximately equal to the GZK-scale l,13l 
of ~ 6 X lO'^eV. Agreement in the spectral characteristics between the two cosmic rays incident in the six degree elevation 
range between payload horizontal 6*^=0 and the Earth's horizon at 6^ - -6°, compared to the 14 cosmic rays observed via their 
reflections, indicate that surface roughness effects are not large, albeit with very limited statistics. That sample was supplemented 
with an additional four reflected events observed with ANITA-II. Interestingly, two of the ANITA-II events featured reflection 
points consistent with sea ice, for which the dielectric contrast relative to the air, and therefore the fractional reflected power, is 
expected to be larger compared to reflections from sheet ice. Using the coordinate convention to be followed later in this paper for 
comparing intensity measurements to expectations from the Fresnel equations, with the signal angle of incidence 0, defined relative 
to a tangent on the surface, the incidence angles for the ANITA radio-reflected EAS sample is shown in Figure |2] From simple 
solid angle considerations {dQ. ~ cos OdO), we expect the acceptance to decrease with elevation angle. Additionally, since the 
antennas are canted at -10 degrees relative to the payload horizontal, large values of 6 are increasingly experimentally suppressed. 
Correspondingly, the observed EAS sample is primarily detected in the angular range 6 <25 degrees. 
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Figure I : Interferometric map for simulated 200, 400, 
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Figure 2: Elevation angle 6, , relative to the horizon, for the 14 ANITA-I, and 4 ANITA-II cosmic ray events observed via their radio emissions. Indicated distances 
are estimated. 
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107 2. Surface Roughness, General Considerations 

108 According to the Rayleigh criterion, the magnitude of variation in surface height h, as measured along the Hne of the incoming 

109 radar must satisfy the condition hcos(<p) < 0.125 (0/, here denotes the "look angle") for the surface to be considered 'smooth'. 

110 In the limit of perfect smoothness (h — » 0, corresponding to "specular reflection"), the Fresnel equations prescribe the fractional 

111 reflected signal intensity, as a function of incidence angle Oj. Deviations from smoothness result in a ^/-dependent redistribution 

112 of the reflected signal power. Depending on the geometry of the surface features, the wavelength, and the dielectric constrast 

113 across a surface, surface roughness can translate into either an enhancement, or reduction in the magnitude of reflected power, and 

114 correspondingly an anti-correlated reduction or enhancement in the amount of power transmitted across the surface. 

115 3. Antarctic Surface Reflectivity 

116 3.1. Visible Frequencies 

117 The polar albedo is one of the crucial inputs to climate models, given its importance in determining the thermal equilibrium 

118 between Earth and Sun. Since the Solar power spectrum peaks in the visible, measurements of Solar reflectivity have typically 

119 focused on the 3900-6600 A regime| 14|. At such scales, sastrugi effects dominate, as a function of look angle (pi relative to the 

120 sastrugi alignment. Long-term studies (particularly at South Pole Station) have established periodic variations in the Solar albedo 

121 in the visible wavelength interval, over ~5 year time scales| 15 1. Such measurements can also be used as a proxy for wind direction 

122 measurements lfT6l . to the extent that albedo observations probe wind-blown sastrugi orientation. 

123 3.2. Prior Measurements at Lower Frequencies 

124 The Earth has, to date, been extensively mapped in the Ku-band (12-18 GHz) lfT7l [TSll and also in the lower-frequency C- 

125 band {2-A GHz), most recently, and most comprehensively by the Canadian-based RADARSAT satellites. These relatively high 

126 frequencies, and correspondingly small wavelengths compared to the typical scale of surface features probe not only pure 'surface' 

127 features, but also sub-surface volumetric scattering due to, e.g., buried sastrugi. The absolute(/relative) RADARSAT signal return 

128 strength calibration of 2 dB(/l dB) is sufficient to map the backscattering power across the Antarctic continent (Fig. [3]l, which 

129 exhibits variations of order 20 dB. Those results are qualitatively consistent with the Envisat satellite Ku-band backscatter[T91, 

130 albeit with some large deviations, particularly around the Ross Ice Shelf. 

131 The Envisat C-band backscattering results generally follow the trend shown in Figure |4] but with the C-band backscattering 

132 strength over West Antarctica reduced by ~2-5 dB; the backscattering signal in East Antarctica is in general agreement with the 
Ku-band scattering to within 2-3 dB. 



backscatter in ku-band from Envisat altimeter 




Maximum contrast corresponds to a variation of almost 20 dB in 
measured reflected power] 17 1. No corrections have been made for 
look angle eft'ects. 

Figure 4: Evisat Ku-band backscattering coefficientfl91, in units of 
dB. No corrections have been made for look angle effects. 

133 

134 Comparable data are not readily available at the higher frequencies of interest to the ANITA mission, although the precision 

135 of three-dimensional satellite maps of the entire Antarctic continent, such as those provided by LandSat, is now approaching 

136 dekameter scales. Proposals are now under consideration to convert photographic images collected during traverses of the Antarctic 

137 interior into surface roughness maps, although those would also sample only a very limited fraction of the Antarctic terrain. 

138 Using the dedicated Taylor Dome discone transmitter, a previous ANITA-I study found reasonable agreement between the 

139 measured radio frequency signal strength at the gondola and the expectation from Snell's Law, knowing the relative geometry of 

140 ANITA-I and the Taylor Dome calibration pulser (Figure |5]l, albeit over a very limited angular range. That study indicates that 
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Figure 5: Measured ANITA-I calibration pulsar received signal voltage (data points), compared to signal voltages without inclusion of Fresnel factors (labeled 
"a/r" and "Ao/r"), and with inclusion of estimated Fresnel transmission coefficient across ice-air interface ((Ao/r)*Fresnel). 



41 surface roughness effects within that range, and in that portion of the continent, are not large. A direct measurement of surface 

42 height variation[i2| at the Taylor Dome calibration pulser site, as well as studies of local radio-frequency propagation through the 

43 surface|3 |, also implied only modest surface roughness RF effects, for that one geographically limited site. 



44 4. Interferometry Applied to Solar Imaging 

45 Absent bursts or flares, the Sun conveniently provides the ANITA mission an uninterrupted, unpolarized, (relatively) constant 

46 power calibration source. However, since the Solar emissions in the ANITA bandpass are typically sub-threshold relative to the 

47 ANITA trigger requirements, the Solar image is only clearly visible after adding many interferometric images in an heliocentric 

48 coordinate system. For our current analysis, we add ~ 10^-10'* events, excluding triggers containing high signal-to-noise sources 

49 to elucidate the Solar signal. 

50 To illustrate the expected signals. Figures |6] and |7] show simulated interferometric maps, for which a single point source, with 

51 a flat frequency spectrum, has been modeled at the indicated elevation angle. Although reflection effects are absent from these 

52 maps, we nevertheless observe some fringing characteristic of ANITA interferograms, and resulting from regular lateral or vertical 

53 spacing of the ANITA horn receiver antennas. We note that for a source at a pay load elevation 6^, fringing tends to produce 

54 an enhancement at an elevation angle typically a few degrees below -6^ and therefore contributes some contamination to our 

55 measured reflection power, with a magnitude that grows with 6^, requiring a small correction to our raw, extracted reflection 

56 signals. 

57 Fig. [8] shows the image of the Sun both in vertical and horizontal polarizations in data, after co-adding 10000 events, acquired 

58 over a time of approximately 20 minutes. During that time, the Sun typically moves 0.5 degrees in the sky. As the Sun has no 

59 intrinsic polarization, we expect approximately equal source intensities in each polarization, consistent with observation. We note 

60 the presence of fringe effects in the interferometric map; such background effects are additionally complicated by the celestial 

61 proxmity of the galactic center to the Sun itself. Over the course of the ANITA-II flight, Sagittarius A* is approximately 5 degrees 

62 higher in elevation relative to the Sun, and within 25 degrees in azimuth, crossing the Solar azimuth roughly midway into the 

63 flight. 

64 4.1. Cross-check of measurement procedure 

65 We have conducted cross-checks to verify the numerical veracity of the interferometric technique, as applied to the Sun, as 

66 follows. 



67 4.1.1. Variation of Solar signal strength with time 

68 For a stationary summer observer at latitudes below the Antarctic circle, the Sun is expected to trace a sinusoid across the sky 

69 over a diurnal period. We observe an obvious 24-hour cycle of Solar source intensity (Figure|9|l, as tracked for an 8-day period of 

70 the flight following the balloon's launch from McMurdo Station. 



71 4.7.2. Antenna Gain vs. Elevation 

72 Over the course of a typical one-month Dec. -Jan. Long-Duration Balloon (LDB) Antarctic mission, the Sun's elevation in the 

73 sky, relative to the payload "horizontal" Opayioad, varies between 10 and 40 degrees, and therefore allows an in-flight determination 

74 of the polar beam pattern for the ANITA dual-polarization quad-ridged horn antennas. Figure [T0| overlays the dependence of the 

75 peak interferometric Solar power (defined as the magnitude of the highest-amplitude pixel in the Solar interferometric map) on 

5 
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Figure 6: Interferometric map for simulated source at 10, 15, and 
20 degree payload elevations, and projections onto elevation axis 
(right). Color scale has units of (Voltage)- . 
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Figure 7: Interferometric map for simulated source at 25, 30, and 35 
degree payload elevations (left), and projections onto elevation axis 
(right). Color scale has units of (Voltage)^ . 



VPol 




-150 -100 -50 50 100 150 



HPol 




-150 -100 -50 50 100 150 



Figure 8: Interferograms (top, VPol and bottom, HPol) formed by tracking the location of the Sun for 10000 events. Normalization is identical for both top and 
bottom plots. For this plot, we use payload coordinates, for which the horizontal plane is defined as 0°, with elevation angles towards the zenith/(nadir) defined 
as positive/(negative). Interferograms are stacked to account for both the motion of the Sun, as well as the translational motion of the ANITA gondola across the 
Earth's surface and the rotational motions of ANITA relative to its vertical and horizontal axes. Horizontal fringes more evident in VPol are consistent with effects 
due to the warm/cold ice/air discontinuity at the horizon, indicated by solid black line at 6'' = -6 degrees. 
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Figure 9: Solar signal strength as a function of number of hours into flight. Normalization is identical for both top and bottom plots. 



Solar elevation, compared with the expectation from a simple Gaussian parametrization of the Seavey horn beam pattern. We 
observe that the Solar interferometric map can be used, in-flight, to calibrate the off-axis beam pattern to an accuracy of a few 
degrees. 

4.1.3. Solar Power as a function of Frequency 



Figure 11 shows the absolute spectral power measured at two terrestrial observatories l|20l 1211 . averaged over the time period 
comprised by the ANITA-II mission. Overlaid on that graph is the direct Solar power, as measured interferometrically by ANITA- 
II. We observe reasonable agreement in the shape of the two spectra. 



5. Measurement of Surface Radio Reflectivity K 

5.1. Geometry 

As indicated in Figure |2] the typical displacement between the ANITA balloon and the surface reflection point is of order 



hundreds of km. For a sub-sample of data spaced out over the duration of the ANITA-II flight. Fig. 12 shows the latitude and 
longitude of both the gondola and also the ice surface reflection location on the ice surface calculated using the known gondola and 
Solar positions at the time of each event trigger, and taking into account Earth curvature effects. The ANITA flight path samples 
the interior ice sheet, as well as shelf ice and some sea ice at latitudes lower than the Ross Ice Shelf. 



202 
203 
204 
205 
206 



5.2. Quantitative Technique 

Similar to FigurejSj we can create an enhanced radio frequency reflected Solar image by adding many events. However, whereas 
Figure [8] is obtained by correcting each interferogram, event-by-event, for both the balloon motion as well as the celestial motion 
of the Sun to maintain coherence of the direct Solar image, proper accounting of the reflection requires correcting for balloon 
motion as well as the motion of the Solar reflection point across the ice surface. Although the location of the Solar reflection 
on the surface is obviously closely correlated to the direct Solar motion in azimuth, and also closely anti-correlated to the direct 
Solar motion in elevation, corrections due to local Earth curvature, surface elevation and the viewing angle relative to the surface 
can easily smear the reflection image by an amount exceeding the intrinsic expected Solar beam spot size (~1 degree) and must 
therefore be carefully tracked as the reflection images are stacked. We note that large-scale variations in surface slope gradients in 
Antarctica typically occur over hundreds of km, and constitute, at most, a 0. 1 degree smearing for the analysis presented below. 

For this analysis, to determine the relative reflected:direct observed signal strength, we fit the direct and reflected peaks to 
Gaussian functions, and take the integral under our Gaussian signal parametrization as our primary estimate of the signal power. 



Figures [13] and 14 show the fits to the one-dimensional projections of the reflection-stacked interferometric maps corresponding to 
Fig. [8] from which we extract signal estimates. For simplicity, we have centered the Solar image in the center of the graphs shown. 

As a cross-check, we have also performed a second, independent analysis, which uses the maximum amplitude pixel in the 
interferometric map, in the vicinity of the known Solar, or reflected location in the sky, as a measure of the signal intensity. This 
technique works well for cases of high signal to noise, which is true for the HPol reflected signal, but becomes decreasingly 
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Figure 10: Parametrized antenna gain of varying beam width (lines) 
compared to relative Solar intensities extracted from Hpol (top) 
and Vpol (bottom) inteiferograms. Color scale indicates number 
of events observed at a given point in graph. 
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Figure 1 1 : Solar intensities measured at two Solai' observatories 
in Learmonth, Australia (open circles) and Trieste, Italy (crosses) 
compared to the frequency banded Solar intensities retrieved from 
the interferograms (Hpol, open triangles and Vpol, open squares). 
The Trieste data are corrected for latitude of the observations. The 
ANITA data are obtained from interferograms filtered into five 
bands (200-400 MHz, 400-600 MHz, 600-800 MHz, 800-1000 
MHz, 1000-1200 MHz); points are placed at the center frequency 
of each band. For each band, we use the same sample of 8000 co- 
added event triggers. Normalizations have been set to be approxi- 
mately equal at 200 MHz. 
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Figure 13: Sample fit to projection of HPol, Sun-centered Solar 
interferometric map onto elevation axis, after cutting ±10 degrees 
around (fi = 0. Fit includes two Gaussian-shaped signals, corre- 
sponding to direct Solar signal and reflected surface reflection, plus 
second-order polynomial to approximate non-zero background un- 
der Gaussian signals. 



Figure 14: Sample fit of projection of VPol, Sun-centered interfer- 
ometric map onto elevation axis, after cutting ±10 degrees around 
= 0. Fit includes two Gaussian-shaped signals, corresponding 
to direct Solar signal and reflected surface reflection, plus second- 
order polynomial to approximate non-zero background under Gaus- 
sian signals. 
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reliabile for low signal-to-noise; correspondingly, we have only applied this cross-check to the HPol case. Fig. 15 displays the 
average ratio K of the Solar radio reflection power relative to the direct Solar image, for both vertical as well as horizontal antenna 
polarizations, including all antenna gain connections, as a function of the reflection angle of the Solar image off the Earth's surface, 
relative to the balloon. Also shown in Fig. [TSjare the expected Fresnel reflection power coefficients, assuming specular reflection 
off the Antarctic surface, and taking the index of refraction of the surface snow to be risurface - 135. We note that the gain 



212 corrections at large incidence angle become quite large for the direct Solar signal, as indicated in Figure 16 



We also note that the vertical spread of points relative to the Fresnel expectation in Figure 15 is approximately a factor of three 



larger than the vertical spread of points shown in Figure 10 indicating that the systematic errors for extracting the signal amplitude 
of the weaker surface reflection are considerably larger than the signal estimation uncertainties inherent in the interferometric maps 
themselves. 

5.3. Interpretation 

Overall, we observe reasonable agreement between the experimental points and the expectation for specular reflection. However, 
we observe an apparent excess, most pronounced for VPol, of the received signal relative to the expectation from the Fresnel 
equations, at higher angles. In particular, we do not observe the predicted signal dimunition (approaching zero signal strength) 
at the Brewster angle, which can be attributed to some combination of a) the result of surface roughness effects, which result in 
non-zero signal at the Brewster angle, b) under-accounting of the effects of fringing, c) the asymmetry in the ANITA-II event 
trigger, which tends to favor thermal excursions of VPol over HPol, and d) an under-accounting of the cross-polarization isolation 



between VPol and HPol. Such cross -polarization effects are, in principle, accounted for in Figure 15 but with the assumption 
that the cross-polarization is constant as a function of frequency and off-boresight angle. Our results suggest that isolation may be 
poorer at large off-boresight angles and increasing with frequency. 

Rigorous simulation of surface effects is beyond the scope of this current work, requiring more extensive data samples on the 
centimeter-scale and shape of the surface features, as well as a cpu-intensive integration across the Antarctic surface. We here 
note only that a simple-minded two-dimensional model comprising an average surface slope of 0.2 radians with a Gaussian power 
spectrum can, at least, reproduce the essential features of our VPol observation, favoring a slight deficit of signal at low incident 



231 angles relative to the basehne expectation, and an excess at larger incidence angles (Fig. 17 1 



232 5.4. Cross-check: Comparison of reflected intensity over sea ice vs. shelf ice 

233 The reflected signal signal strength at an interface depends on the dielectric contrast between the air and the reflecting layer. The 

234 Solar image reflected off sea ice, having relatively higher salt content, should be brighter than the Solar image reflected off' sheet 



ice. As illustrated in Figure 18 we correspondingly note the largest enhancement in the reflected signal strength in the coastal 



region beyond the Ross Ice Shelf, consistent with expectation. 
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Figure 15: Ratio of ANITA-II measured reflected Solar power to directly-received Solar power, as a function of incident elevation angle, for both HPol and VPol 
antennas. Overlaid are also calculated values as obtained from direct application of the Fresnel equations to calculate the expected reflected signal power K. 
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Simulation (green) vs. VPol Fresnel (black) 




Elevation angle 6, wrt surface (deg) 

Figure 16: Magnitude of applied gain correction, as function of ele- 
vation angle, for direct Solar signal as well as reflection signal. 

Figure 17: Simulation of expected VPol reflection signal strength 
using simple model of Gaussian surface of inclination angle 0.2 ra- 
dians, and with random scattering at all points. Simulation is qual- 
itatively consistent with VPol measurement, showing deficit at low 
incidence angles and excess at higher incidence angles. 




Figure 18: Deviation between measured HPol reflection coefficient 
relative to reflection coefficient expected at an air-ice sheet interface, 
as a function of position on the Antarctic continent. 



Figure 19: Deviation between measured VPol reflection coefficient 
relative to reflection coefficient expected at an air-ice sheet interface, 
as a function of position on the Antarctic continent. 
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237 6. Discussion and Summary 



To determine the energy of a typical down-coming air shower event, one must estimate the reflectivity of the Antarctic ice 
surface. The reflectivity is, in turn, closely related to the physical scale of surface features, which depend on such parameters as 
local wind speed and direction and surface elevation gradient. Surface roughness also has implications for neutrino detection, as 
high-frequency, small-wavelength radio waves resulting from in-ice neutrino interactions may decohere as they pass through the 
ice-air interface, en route to detection at the ANITA gondola. The wavelength regime to which ANITA is sensitive (125 cm— »25 
cm) approaches the expected size of surface inhomogeneities, and therefore underscores the importance of direct surface effect 
measurements. 

In our current analysis, the ANITA-II balloon-borne radio interferometer has been used to measure the Antarctic surface re- 
flectivity over the ANITA passband. We find reasonable agreement, for both Vertical as well as Horizontally polarized signal 
power, with expectations based on direct application of the Fresnel equations. This analysis therefore provides some basis for 
directly applying the Fresnel coefficients to infer the energies of cosmic rays measured by their surface-reflected radio-frequency 
emissions, as well as the energies of neutrinos detected via their emissions across the air-ice interface. Ongoing analysis will 
incorporate these results, combined with improved simulations of the radio frequency EAS signal, to give revised estimates of the 
observed ANITA-I cosmic ray events. As an indication of what the associated systematic errors in those energy estimates are. 
Figure |20]presents the uncertainty implied by the current analysis. We observe that the typical event-by-event energy uncertainties 
due to surface reflectivity is of order 30-40%, and increasing as the magnitude of the antenna gain correction increases, while the 
deviation between the average reflectivity and expectation, for the large ensemble of events expected for ANITA-3, is typically a 
factor of three smaller 

We stress, however, that for the much higher signal-to-noise typical of EAS events, signal energy estimates can be made directly 
from the event waveforms, without relying directly on the interferometric maps to infer signal strength. Variations in the event-by- 
event energy estimate should be considerably smaller than those encountered in the Solar reflection analysis; the systematic errors 



implied by Figure 20 therefore represent a (very) conservative upper limit. As illustration. Figure 21 shows the signal-to-noise 
power ratio measured for higher amplitude calibration pulser events, indicating an intrinsic event-by-event energy uncertainty a 
factor of ~2 smaller than what our analysis, based on the lower amplitude Solar signal, might otherwise indicate. 



Energy Systematic Error 



SNR for Taylor Dome Calibration Pulser Events 
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Figure 20: Implied event-by-event radio air sliower energy uncer- 
tainty, assuming a) errors are equivalent to tliose obtained in current 
analysis, and b) energy estimate is dominated by HPol response. 
Values are dominated by spread in values in Figure [TI] rather than 
observed deviation between measurement and expectation. 
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Figure 21: Signal-to-noise ratio measured for calibration pulser 
events using in-ice transmitter discone at Taylor Dome, as a function 
of ANITA-II event trigger number. As ANITA-II moves, optimal 
signal reception occurs approximately at event 1,600,000. Typical 
scatter within each bin is observed to be of order 20%, taken to be 
representative of the intrinsic interferometric energy resolution un- 
certainty for high SNR events. 
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